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Abstract 

The origin of continuous energy spectra in large disordered interacting 
quantum systems is one of the Itey unsolved problems in quantum physics. 
While small quantum systems with discrete energy levels are noiseless and 
stay coherent forever in the absence of any coupling to external world, 
most large-scale quantum systems are able to produce thermal bath and 
excitation decay. This intrinsic decoherence is manifested by a broadening 
of energy levels which aquire a finite width. The important question is 
what is the driving force and the mechanism of transition(s) between two 
different types of many-body systems - with and without intrinsic deco- 
herence? Here we address this question via the numerical study of energy 
level statistics of a system of spins- i with anisotropic exchange interactions 
and random transverse fields. Our results present the first evidence for a 
well-defined quantum phase transition between domains of discrete and 
continous many-body spectra in a class of random spin models. Because 
this model also describes the physics of the superconductor-insulator tran- 
sition in disordered superconductors like InO and similar materials, our 
results imply the appearance of novel insulating phases in the vicinity of 
this transition. 

Although quantum phase transitions between regimes with and without in- 
trinsic decoherence were first discussed by Anderson in 1958 [5] , their under- 
standing is still incomplete. More recently, a number of works have studied these 
transitions in strongly disordered electron systems with weak repulsion [Tl [T71 H] 
and reached the conclusion [T7l|4] that a finite-temperature phase transition takes 
place between two regimes, "weakly insulating" and "strongly insulating", that 
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are characterized by non-zero and zero conductivities respectively. In terms of 
the many-particle excitation spectra this result translates into the existence of 
an extensive energy threshold, with a critical excitation energy £c oc VT^ that 
separates the discrete spectrum of excitations at E < £c and continous one at 
E > £c for any large but finite subsystem with V degrees of freedom. It also 
implies that the low-temperature state is free from decoherence because the life- 
time of all excitations is infinite. The phase transition predicted in the work [4] 
was not observed yet; instead, in hopping insulators one usually observes that 
conductivity vanishes continously at T — >■ 0: a(T) ^ q-{t' /Ty ^ with a < 1. 

Different physical properties of the excitations at low and high energies in the 
infinite system are reflected in different statististical properties of the spectra 
of finite systems sX E < £c and Sit E > £c- Intuitively, if the eigenvectors 
are extended, as expected for the state where local excitations decay, they are 
subject to inter-level repulsion. Conversely, if the eigenvectors are localized, 
eigenvalues corresponding to excitations localized in different parts of the system 
are independent and one expects a Poisson distribution of energy levels. 

Let us give a general qualitative argument in supporting this statement. 
Consider a small perturbation of the Hamiltonian that controls the dynamics of 
a generic quantum system in thermodynamic limit: 

H^H{l + 4>{t,x)) (1) 

where <j){t, x) is a generic slow function of coordinates and time. A small pertur- 
bation of this type results in the slow (adiabatic) motion of energy levels En{t). 
In the absence of level repulsion, different levels cross without affecting each 
other, so that this motion leads only to the total phase of the wave function. 
Because the field (j){t,x) (which is similar to the gravitational potential [19]) is 
conjugated to the energy density, the absence of response to it implies absence 
of the energy fiux. An excitation with energy AE localized around point x 
acquires phase exp{—iAE<j){t,x)) due to perturbation ([T]); in contrast, a de- 
localized excitation becomes a superposition of other excitations. Thus, the 
absence of the response also implies that excitations are localized and do not 
decay. We conclude that the absence of level repulsion implies the localization 
of excitations, absence of their decay and of the energy fiux, i.e. formation of 
a strong insulator. Because level statistics can be studied for relatively small 
systems, this correspondence between level statistics and physical properties in 
thermodynic limit provide a convenient numerical tool to predict the properties 
of physical systems. This strategy has been used for instance in |20j . 

One of the best experimental systems to test these general ideas is pro- 
vided by strongly disordered superconductors (InO, TiN) in the vicinity of 
the disorder- induced superconductor- insulator transition [28] ; for recent reviews 
see [16] and the introduction of the paper [11] . One advantage of these systems 
is that the transition can be fine-tuned by a magnetic field. In the insulat- 
ing state they demonstrate purely activated behaviour cr{T) ^ e~^'^ with 
a = l.[371[3T] Far in the insulating state [5H] this behaviour can be understood 
in terms of a single electron pseudogap [351 HI] that results from binding of 
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localized electrons [121 [H]- However, the persistance of an activated behavior in 
the vicinity of the superconductor-insulator transition (SIT) where one expects 
the appearance of low energy Cooper pairs requires another explanation. This 
was proposed in the work [TSl [13] in terms of a collective mobility edge Cc which 
separates the domain of localized excitations with energies w < ec from the 
domain of delocalized modes with w > Ec, which serves as an intrinsic thermal 
bath. The threshold ec{g) depends on the parameter characterizing disorder, g. 
It vanishes at g = gc at which superconductivity appears. At g < gc the thresh- 
old energy Cc is intensive [HI [13] (it does not grow with the size of the system). 
In this regime, the density of excitations is proportional to exp(— ec/T), leading 
to the activated behaviour of conductivity. At even stronger disorder, g = g*, 
the threshold energy e^g) diverges. 

Logically there are two possibilities; the first one is that the divergence of the 
threshold ec at — )• g* is cut off by the system size, so that at stronger disorder 
Cc scales with the system volume, V. In this situation one expects the finite 
temperature transition between insulator and hard insulator predicted in |17[ H]. 
In the second scenario, the energy scale Cc becomes really infinite at 5 — > g*, 
which implies a complete localization and the absence of a thermal equilibrium 
at all temperatures. We will show below that both scenarios can be realized 
in discrete spin models depending on the detailed form of the Hamiltonian: in 
the models with purely local density-density interaction Cc is infinite at g < g*, 
but the addition of non-local density-density interaction terms results in the 
appearance of a narrow range gi < g < g* where Cc is extensive, signaling the 
finite temperature transition in this regime. This result raises the possibility 
that a finite temperature transition between strong and weak insulators might be 
observed in InO, TiN or similar films very close to the superconductor- insulator 
transition. 

The theoretical work [TH] [13] is based on several approximations. First, it 
assumes the presence of a strong pseudogap, which allows one to reduce the 
Hilbert space of the full electron problem to the one spanned by "pseudospin"-i 
variables s^,Sj, that describe [3] creation-annihilation and counting operators 
of the localized Cooper pairs. This assumption is borne out by experimental 
results [25]. On the theoretical side, its origin lies in the fractal stucture of the 
wave functions near to the mobility edge of the single electron problem [l2l [TTj . It 
is further supported by direct numerical simulations of electrons with attractive 
interactions in a random potential [7]. The second approximation is the use of 
a recursion-equation technique, which becomes exact only for a special tree-like 
structure of a lattice completely devoid of small loops (Bethe lattice). Finally, 
the formalism developed in these works does not include explicitely collective 
excitations of high energy involving many spins. Therefore this formalism does 
not allow to completely settle whether the mobility edge is intensive or extensive. 
Note that although another prediction of our formalism, namely the anomalous 
broadening of the superconducting order parameter distribution in the vicinity 
of transition, was confirmed experimentally [25] , no such data exist for insulating 
phase. 
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In the present Letter we provide a direct numerical proof of the validity 
of the many-body localization scenario developed in[T8l [13] for a more realistic 
random spin lattice and taking into account the whole excitation spectrum. Our 
basic model contains spins ^ that describe the absence or presence of Cooper 
pair on a localized single electron state. These spins are subject to random 
fields along the z direction and are coupled by XY interactions The main 
physical result of our study are the phase diagrams, shown in Figs. [T]and[2l in 
terms of the dimensionless transverse spin coupling g and excitation energy w 
(Fig[T]) and temperature (Figl2]). A size-independent (non-extensive) threshold 
energy ec{g) is found in this "minimal" model, Eq.(l2]), for g € {g*,gc), while 
all eigenstates are localized a.t g < g*. This implies that, for this Hamiltonian, 
the transition line separating the weak and hard insulators does not depend 
on temperature. On the other hand, we shall show that, in a model where 
the Hamiltonian is modified and includes additional interactions between the z 
components of spins, the threshold energy remains finite but extensive in finite 
systems in some range of 5 < g* , see Fig. (3] This leads to the temperature 
driven transition similar to the one predicted in [T71 0] . 

The phase diagrams shown in Figs. [T]and[2]have another important feature: 
they predict a direct transition between superconductor and insulator with char- 
acteristic energy scales that go down to zero continously on both sides of the 
transition, similar (but different in the important details) to the transition ex- 
pected in dual theories [TH [TS]. In the superconducting state this energy scale 
is given by the typical value of the order parameter or the transition tempera- 
ture, while in the insulating states it is the value of the threshold energy ec{g), 
which implies the Arrhenius behavior of the resistivity at very low temperatures. 
These predictions are in agreement with the detailed studies of the transition 
driven by a magnetic field in films [25] and in Josephson junction arrays with 
moderate Ej /E^^^. 

As discussed in details elsewhere [TTl fT3] a superconductor with a large pseu- 
dogap and a weak long-range Coulomb repulsion is faithfully described by the 
model of spins- ^ : 

H = -Wj2^.s^, - J2 + sTsD (2) 

' (ij) 

where points i,j belong to a random graph G with a fixed coordination number 
Si — ifji are spin-i operators, sf = sf ± is\, the sum goes over all 

different pairs of nearest neighbors i,j on G, and all nonzero matrix elements 
are equal to J^^ = Wg/{Z — 1). The random energies are uncorrelated 
at different sites and choosen from the box distribution Q(^) = — |^|), 

corresponding to bandwidth W . 

Eigenstates of the Hamiltonian ([2]) are vectors in the 2^-dimensional Fock 
space (N is the total number of sites of G). As discussed above, the transition 
between decoherent and coherent states can be deduced from the change in the 
statistics of the exact eigen-energies Ei of the Hamitlonian. To identify the 
energy level statistics, we study the dimensionless parameter r„ G [0, 1] defined 
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Fig. 1 : Phase diagram of a strongly disordered superconductor with a large pseu- 
dogap as a function of interaction constant g. The full lines show the 
predictions of the analytical theory of the model ([2]) for the critical tem- 
perature (right vertical axis) and the threshold energy, e, (left axis) of 
spin flip excitations in infinite random graph with Z = 3 neighbors. The 
vertical ovals show the values of the critical coupling constant that corre- 
spond to a transition between different types of spectra for different ener- 
gies E in the finite random graph model of a small size {N = 16 — 20) as 
determined by direct numerical simulations. The uppermost oval shows 
the transition at the many-body band center (corresponding to E 1) 
that sets a lower bound for the critical g{E). The thick dashed line 
shows the position of the spectral threshold for single-spin excitations 
with energy e adjusted by finite-size effects, as explained in the main 
text and in the Supplement IA.21 The small circles show the typical en- 
ergy of the single-spin excitations, e{E), that give the main contribution 
to the many body excitations studied in direct numerical simulations. 
The good agreement between their position and expectations (dashed 
line) confirms the validity of the cavity method [THl [13] that is used to 
obtain the results in large systems. The very small change in the criti- 
cal value of the coupling constant between excitations at energy E « 2.0 
and the band center implies that all excitations, at high and low energies, 
become localized when g < g* ■ 
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Fig. 2: Phase diagram in the temperature - coupHng constant plane for the 
model (jl} with Z — 3 {K = 2), obtained from the solution of cavity 
equations. The strength of the s^s^ interaction is — 0.1. This in- 
teraction leads to a qualitatively new phenomenon, the appearance of a 
finite temperature transition between weak and strong insulators. In the 
weak insulator, excitations at sufficiently high energies can decay even 
at zero temperature. A non-zero temperature results in non-zero relax- 
ation of all excitations, even the ones of lowest energy. In contrast, in 
the strong insulator, no excitation with intensive energy can decay. As 
the interaction constant is decreased, the temperature separating these 
phases goes to infinity at g = gi. At smaller coupling g < gi, all exci- 
tations, even those with extensive energy remain localized. The value of 
gi « 0.042 is approximately equal to 0.30gc- The ratio gi/gc = 0.3 is in 
good agreement with the results of the direct diagonalization on small 
graphs, as can be seen from the critical values of J^^ found in FiglSl 
0.02 for the middle of the band band and 0.74 for the low energies. 
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Fig. 3: The energy- level statistics is characterized by the average (r„) that dis- 
tinguishes Wigner-Dyson and Poisson distributions (values of (r^) cor- 
responding to these distributions are shown by dashed lines). The left 
panel shows the statistics of the low-energy excitations in the energy in- 
terval {Egsi ^-gs + l.S) as a function of the transverse interaction constant, 
J^**, for the Z = Z random graph with bandwidth W — 1. The mid- 
dle panel shows similar results for intermediate energies, and the right 
panel corresponds to high energies, close to the center of the many-body 
spectrum. 



as 

r„ = min((5„, dn-i)/ max((5„, 6n-i) (3) 

where Sn — En — En-i and En is the n-th energy level. The average value 
r = (r„) is equal to 0.38 for Poisson level statistics, and to 0.53 for Wigner- 
Dyson statistics. In the limit of infinite systems one expects a sharp transition 
between these two values as a function of g, or of the energy. In a finite system 
the parameter r increases smoothly, but the curve r{g) is expected to become 
steeper as the size is increased. 

The appearance of a unique crossing point of these curves for different sizes 
N implies a well-defined phase transition in the limit of infinite-size systems. 
This is seen in the exact diagonalization of the Hamiltonian ([2|) with Z = 3, 
shown in Fig. [31 Apart from the persistance of the transition for small system 
sizes, these data also show that the critical value of J^^ changes very little as 
the energy is increased to the center of the band. This proves that high energy 
states in this model become localized together with the low energy ones, in more 
physical terms it implies that no new decay channels appear at high energies. 
The absence of new decay processes at high energies is the physical reason why 
in this model one does not get the intermediate phase with extensive energy 
threshold. Also, these data show that the transition happens at a value rn{g) 
that is close to its value 0.38 expected for Poisson statistics. This is due to a 
large distribution of relaxation rates in these systems [13] that implies that in 
case of small systems many realizations of the random energies give localized 
states, while delocalization only happens with a small probability as we explain 
in a more detail below. 

In order to compare quantitatively the results of the direct diagonalization 
shown in Fig. [3] with the predictions of the theory '13 we need to take into 
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account the finite size effects. They are very significant for the sizes (number of 
spins N < 20) available for direct diagonalization, for two reasons. First, in a 
finite system the crossover from Poisson to Wigner-Dyson statistics is expected 
to occur when the level spacing S{E) = h'~^{E) becomes approximately equal 
to the level width r{E) expected theoretically. The latter quickly becomes ex- 
ponentially small [T5j when the transition in the infinite system is approached, 
so that the condition S{E) = r{E) in a finite (not very large) system is satisfied 
far away from the transition in the infinite system. Second, the level widths 
fluctuate wildly from one finite system to another in small systems. Because 
the infinite system can be viewed as being composed of small ones, each of these 
parts having many neighbors, even a small probability to find a delocalized (fi- 
nite width) level in a small system is sufficient for delocalization in the infinite 
system. This makes the typical decay rate in the infinite system much larger 
than the one in a collection of small ones. Both effects combine to push up the 
apparent critical value of g by roughly a factor of two. Finally, for the quanti- 
tative comparison, we need to take into accout that the critical point condition 
6{E) = r{E) applies to the level spacing and total level width in the many body 
systems, whereas the analytical theory |13) gives the level width of individual 
spin fiips. The typical level at energies E > W consists of a few spin fiips, 
which decay rates add to the total width r(£'). The main contribution to the 
typical level at energy E >W comes from single-spin excitations in a relatively 
narrow energy window e{E) « ^ EW/N, see Supplement lA.il This allows us to 
interpret the crossing point in exact diagonalization as the delocalization of in- 
dividual spin fiips as shown in Fig. [TJ We present the details in the Supplement 

These results are consistent with the recent findings of paper [8] which stud- 
ied a model similar to our ([2]) but with identical s^s~ -couplings between all 
pairs of spins, allowing exact integrability of the Hamiltonian. The conclusion 
reached in this work is that many-body delocalized state disappears in exactly 
integrable case but reappears when variations of couplings between spins are 
allowed that destroys integrability. In the latter case one expects a sharp tran- 
sition between localized and delocalized regimes similar to the one found in 
model (O. 

We now discuss the generality of the conclusions reached above. The absence 
of an intermediate phase with extensive energy threshold can be traced back to 
the irrelevance of interactions between individual spin fiips. One can thus expect 
new physics to appear if such interaction is introduced. The simplest and most 
physical model that has an additional interaction between spin flips includes an 
additional "longitudinal" spin-spin interaction: 

* (y) 

that can be viewed as due to the long range part of Coulomb interaction ( < 0) 
or to phonon mediated attraction between electrons {{J^ > 0) [TT]. As will be 
clear below, the results do not depend on the sign of ^ W. The first two 
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terms in the Hamiltonian Q describe a classical Ising magnet in a random 
field. We shall be interested in its disordered phase which is realized when the 
longitudinal interaction is small, <C W. In this case the classical eigenstates 
of this magnet coincide with that of independent spins, with energies that are 
weakly modified by the interaction. In the absence of transverse interactions 
the excitations of this magnet are individual spin flips with energies = W£,i + 
^j- the ground state the direction of almost all spins is determined 
by the sign of the random field, so the presence of a small affects weakly the 
density of states of these low-energy excitations and their decay due to transverse 
interaction. Thus we expect that the low-energy properties of the spectrum 
remain similar to the model with zero J^. The situation is very different at 
high energies when the decay of a given spin occurs against the background 
of different spin configurations. Physically, a spin excitation at energy e has a 
much larger chance to propagate through a given site if the energy of the spin 
flip at this site is close to e. The energy of a spin flip is distributed with the 
probability density where ^ = e/W^ if all surrounding spins are locked 

into a fixed configuration. In contrast, at extensive energies, the surrounding 
spins acquire many different configurations. This increases the probability that 
a given site is in resonance with the excitation. This effect can be also described 
as being due to the interaction between individual spin flip processes that allow 
new channels for the decay of these excitations. This should result in the large 
suppression of the critical value of the interaction constant at high energies. 

The exact diagonalization of the Hamiltonian on a Z — 3 random graph 
with — 0.1 confirms these expectations (Figure S]). In presence J^, the value 
(7c(0) is shifted slightly downwards for low energies ((Jc = 0.10 — ^ 0.074). The 
shift becomes substantial for medium energies (Jc = 0.066 — ?> 0.04) and very 
large for energies close to the band center: Jc — 0.06 — ^ 0.02. A slight shift 
of the critical value of the interaction constant at low and medium energies is 
probably due to the fact that at these energies a typical excited state contains 
more than one spin-flip excitation. As explained above, the interaction between 
these excitations leads to delocalization. 

The solution of cavity equations, similar to those solved in [T3] (see Sup- 
plement IA.3|) confirms the appearance of the extensive threshold for ^ 
and gives the dependence of the transition temperature Ti{g) on the coupling 
constant g that we show in Fig. [2] for = 0.1. At g < g* and T < Ti{g) 
all excitations with intensive energies are localized, no transport of any sort is 
possible. In contrast, at higher temperatures T > Ti{g) all excitations, even 
those with low energies, acquire a non-zero width. At g > g* and zero tempera- 
ture, the excitations with low energy are localized while high energy excitations 
decay. This latter distinction is smeared at any non-zero temperature because 
the presence of even a small density of mobile high energy excitations leads to 
a slow decay of even lowest energy ones. 

The main conclusion of our work is the existence of two different insulating 
phases and a prediction of the finite temperature transition between a weak 
insulator characterized by an activated transport and a strong one in which no 
transport occurs. Experimentally, the second transition might be detected by 
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Fig. 4: Level statistics, characterized by the parameter (r„) in the presence of 
weak longitudinal spin coupling J for low energy levels (left panel), in- 
termediate energies (middle panel) and center of the many body band 
(right panel). Even a small coupling = 0.1 has a significant effect, it 
shifts the transition to much smaller values of the transverse coupling g. 



the onset of an anomalously slow thermalization:9j. These results are obtained 
within a lattice spin-i model which serves as a good approximation to the 
description of the superconductor-insulator transition in a number of materials. 
We believe that it might be possible to observe this transition in insulating 
films in a close proximity to the superconducting transition. Due to a finite 
value of the pseudogap in a realistic systems, a temperature-driven transition 
between weak and strong insulating phases will be seen as a sharp crossover in 
resistivity curves R{T), with a rapid growth of the apparent activation energy 
d\nR[T) / d{\/T) as T decreases. Some preliminary experimental evidence for 
this behaviour in TiN was reported in [5], it was also observed in InO^; |23j . 

EC and MF thank the FEDER and Spanish DGI for financial support 
through project FIS 2010-16430. MF acknowledges support through RFBR 
grant 10-02-00554 and RAS program "Quantum physics of condensed matter", 
LI was supported by ARO W911NF-09-1-0395, DARPA HROOll-09-1- 0009 and 
NIRT ECS-0608842. 
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A.1 Density of states 

For the quantitative analysis of the numerical results wc need the energy depen- 
dence of the total density of states v{E) of the Hamiltonian ([2]). As we explain 
in detail below, at low energies i? < it is given by the formula first found by 
H. Bethe [5] for heavy nuclei: 

v{E) = CNeyii>iaVEN) (5) 

where C ^ I and a ^ 1 are functions which depend only weekly on g. This 
equation and the discussion below assume that the band width is chosen as 
W = 1. The numerical data presented in Fig. [5] show that the density of many 
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Fig. 5: Density of states of the Hamiltonian ([2]) obtained for g — 0.04, g = 
0.16 and g = 0.30, as indicated on the panels. The coefficient a{g) is 
determined by the slope of the straight lines that provide the best fit to 
the data points. The upper right panel shows the a{g) dependence. It 
diplays a maximum for g close to the critical point gc- 



body levels of Hamiltonian ([2]) are indeed given by Eq.© for all values of with 
a weakly g-dependent coefhcient a{g). Because the functional form of i'{E) is 
(^-independent, one can understand its origin by considering the simple problem 
of non- interacting spins, defined by the first term of ([2]). 

In this approximation the excitation energy is E = X]i=i n Vi'^i where cr^ 
are equal to or 1, and rji e [0, 1] are quenched random parameters. Let us 
denote by v{E) the typical density of states. Its Laplace transform , for one 
given instance (i.e. one given realization of 77^) is 

D{s) = exp(- 5^ ,7.^.5) = n (1 + ) ■ (6) 

sie{o,i} i 

For a large system, this density of states self-averages, in the sense that log v 



A Supplementary online material 



12 



goes to a well defined limit: 

\ogi>is)=N I dr]\og{l + e-'") . (7) 







In order to reconstruct the typical density of states using the inverse Laplace 
transform, we must find the maximum over s ofEs + N /g d?7log(l + e-^'').The 
value of s at which this expression is maximal satisfies the equation 



N Jo ' 1 + e-'*'' 
leading to 

i^{E) « exp I^V^ ^ dT] In [l + e"^"''] + Esj (9) 

In the limit of small E/N the maximum is found at large s, so that 
E I f°° , a;e'^ c 

— ~ ^ / dx = (10) 

where c = j^- Altogether this shows that the density of states grows as e"^^=-^, 
where 

a:=2VS=^wl.8 (11) 
v3 

For realistic N, however, the minimum is somewhat different as can be seen 
from the Figure |6] in which we show the full expression ([9]) for = 18 and its 
fit to the square root dependence. 

In comparison, the direct numerical simulation of 100 realizations for A = 18 
system shows a similar behavior with slightly smaller coefficient a = 1.25— 1.35. 
The asymptotic ([5]) does not change if we restrict the states to the sector = 
computed numerically, see Figure |6] where we show the results restricted to this 
sector. 

In the energy range 1 < < A^, the most probable number n{E) of single- 
spin excitations whose energies sum up to E, is n{E) « V EN. Typical single- 
spin excitation energies are e{E) = E/N. 

At higher energies, density of states reaches it maximum a.t E = Eo ^ N, 
with lnz/(£'o) ^ A^i where Eq corresponds to the center of the many-body band. 
One can relate extensive energy E = eN of the spin system and its temperature 
T, assuming internal equilibrium: 

T=f^^y' = V) for E«N (12) 
\ dE J a 

Equation (jl2p shows that typical single-spin excitation energy grows with tem- 
perature as e = aT/2. Note that T{E) diverges as energy approaches band 
center, T-^{Eq) = 0. 
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Fig. 6: Density of states and its approximations. The left panel shows the den- 
sity of states given by the full analytical formula ^ and its approxi- 
mation by the simplified square root dependence In = aV NE + b for 
= 18 with a = 1.4 and b = 1.5. The right panel shows the direct 
numerical simulation of 100 realization for the system of 18 spins in sec- 
tor Sz = and its fit to a In = aVnE + b dependence, with a = 1.25, 
b = 0.5. 

A.2 Finite-size effect upon the transition line 

In a finite system the crossover from the Poisson to Wigner-Dyson statistics is 
expected to occur when the states become delocalized. Delocalization implies 
that the width of the level determined self-consistently in the cavity approach 
[13] is of the order of the level spacing. A number of finite size effects make 
the direct comparison between the analytical results for infinite systems [TS' and 
numerical results reported here non-trivial. First, the width, T(E), of the level is 
expected to become exponentially small as the energy threshold is approached: 

r(e) e"--"=(9) 

Furthermore, the coefficient, wi in this equation is much larger (see Ref. |13) ) 
than ec, so the decay rate becomes very small even relatively far from the tran- 
sition line. In a finite system with a not-so-small level spacing, this leads to 
a very significant shift of the apparent critical energy due to finite size effects. 
Second, the width of the levels fluctuate strongly from one graph realization to 
another, see insert in Fig. [T] The large peak at small values of F shows that in 
many realizations the level width is essentially zero. The change in the statistics 
observed in small systems is due to a relatively small number of graphs with 
significant F. This explains why the crossover is always observed at a value of 
r„ that is close to the one of Poisson statistics. In order to compute the criti- 
cal value of the coupling constant expected in finite systems we evaluated the 
probability, P(Fo) to find the level width F > Fq obtained by the solution of 
cavity equations for the systems of this size. The result is shown in Fig. [7|for a 
typical value of e and Fq. The crossing point observed by direct numerical diag- 
onalization is expected to happen when the probability P(Fo = S{E)) becomes 
non-zero. To avoid the problem with numerical errors, we used the condition 
P{To = S{E)) > So with Sq = 0.004 - 0.02 and checked that the results are not 
sensitive to the specific value of So in this interval. Using this condition we take 
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Fig. 7: The solution of cavity equations in finite random graphs gives a widely 
distributed width F of single spin flip excitations. The insert shows the 
probability distribution of r(cL>) for oj = 0.3514^ in a random graph with 
N = 20 spins and a coupling constant g = 0.10. The right-hand peak in 
this distribution corresponds to all graphs in which a significant r(a;) is 
spontaneously formed. The main panel shows the weight of the graphs in 
which a significant T{uj) (defined as InF > InFo) is formed as a function 
of the interaction constant g. The resulting value is only weakly sensitive 
to the value of the cutoff Tq for most w and g. This allows us to determine 
the apparent critical value of the interaction constant g{uj) for the finite 
systems. The results are shown in FiglT]as a dashed line. 




into account both finite size effects and compute the finite-size corrections to 
the infinite-system edg) predicted analytically. The final result for the apparent 
transition for small sizes is shown in Fig. [Hby the green line, it is shifted with 
respect to the infinite-size edg) shown by red line by a factor of two. 

Finally, when comparing the predictions of the cavity equations with the re- 
sults of the exact diagonalization, we have to take into account that high energy 
levels of the whole system correspond to many spin flips: the resulting decay 
rate is the sum of the decay rates of the individual flips. Because of the weak 
(logarithmic) dependence of the probability ^'(Fo) on Fq this fact has a very 
weak effect on the expected transition. Because of the fast (exponential) depen- 
dence of the density of states on the energy, the typical spin flip contributing 
to the high energy level has a well deflned energy given by the saddle point 
solution of section lA.il e{E) ^ y^EW/N. This allows us to map the energies, 
E, studied in the numerical diagonalization to a typical energy of a single flip, 
e{E), and to compare the results of the numerical diagonalization and solution 
of cavity equations as shown in Fig. [1] 
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A.3 Sz - Sz coupling and finite-T phase transition 

In the previous papers |TS1 [T3] we derived a recursion relation for the intrinsic 
level widths of single-spin-flip excitations, assuming a local Bethe lattice 
structure of the lattice: 



( j-y)2 y (13) 



where k(i) are the K ^ Z — 1 neighbours of i in the cavity graph. The same 
recursions can be used for a random graph model studied in the present paper, 
upon neglecting the contributions of closed loops (loops are totally absent on the 
infinite Bethe lattice whereas they are present with parametrically low density 
^ 1/ In TV on a random graph with N sites) . The condition of stationarity of the 
distribution function P{T) generated by the mapping (1131) leads to the solution 
for typical Ttyp{LL>) found in Ref. [13]. 

In presence of sf s| coupling the recursion equations (IT3l) should be modified. 
At nonzero temperature T = it contains the sum over neigboring spin 
configurations with their thermal weights: 

y - — — (14) 

where the internal summation over s^j.^.^ includes all configurations of spin 
variables connected to the spin by z — z links, Zk — e^^"*"'*") are the 

corresponding local partition functions, and = W^|Cfc| + ^n{k) ^n- 

The largest contribution to the sum (jl3l) comes from the site k characterized 
by Cfc that is closest to uj. The same holds for the sum but in addition in this 
sum the value of varies depending on the surrounding spins. This increases 
the probability to find a resonance. To evaluate the importance of this effect 
we consider explicitly the case of Z = 3 (corresponding to the random graphs 
diagonalized in this paper) and the energy at center of the band, uj = W/2. The 
critical value of the coupling constant in this case is determined by the condition 
InS = in the limit of the large system size. Here 



.g\2 1 f e-/3(?,i+«.2) 2cosh/3(^ji -^j-2) e^(«^i+«^-2) 



... 

is the relaxation rate induced deep in the system center by the infinitisimally 
small couplings at the boundary, Zj = 4cosh/3^ji cosh/?^j2, 7 = /W , f3 = 
2T/W and S(,j = — w. Performing the same steps as in Ref. [T3] we average 
over the distribution of ^ using the replica method and we get the condition for 
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the critical value of g{T) 
g — K cxp 



min fix] 

2 X 



X 



g-/3(?i+«2) 2cosh/3(ei -6) e'5«i+«^) 



The inversion of the obtained function g{T) provides the Tj{g) dependence 
shown in Figl2] by the thick violet line. 
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